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Abstract 

We study transverse-momentum-dependent (TMD) gluon distributions and factorization the¬ 
orems for the gluon-initiated semi-inclusive processes at hadron colliders. Gauge-invariant TMD 
gluon distributions are defined, and their relations to the integrated (Feynman) parton distributions 
are explored when the transverse momentum is large. Through explicit calculations, soft-collinear 
factorization is verified at one-loop order for scalar particle production. Summation over large 
double logarithms is made through solving the Collins-Soper equation. We reproduce the known 
result in the limit that the transverse momentum of the scalar is large. 
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I. INTRODUCTION 


In recent years, there has been considerable experimental and theoretical interest in semi- 
inclusive hard processes, from which one hopes to learn important information about the 
nucleon structure and non-perturbative dynamics of quantum chromodynamics (QCD). Rig¬ 
orous theoretical studies in this direction started with the classical work on semi-inclusive 
processes in e+e“ annihilation by Collins and Soper in j^, where a QCD factorization has 
been proved, and nonperturbative transverse-momentum-dependent (TMD) parton distri¬ 
butions and fragmentation functions were introduced Hi The approach has been applied 
to semi-inclusive Drell-Yan production at hadron colliders in and to semi-inclusive deep 
inelastic scattering (SIDIS) in 0,0, in both cases with transverse momentum much larger 
than Aqcd- Meanwhile, spin-dependent TMD parton distributions and fragmentation func¬ 
tions were introduced in various processes, which can generate polarization asymmetries in 
semi-inclusive scattering In the past few years, gauge properties of the TMD 

parton distributions have been investigated 11 MMM, and the result provided a hrm 
theoretical basis for studying the TMD parton distributions. More recently, the factoriza¬ 
tion theorems for the semi-inclusive deep inelastic scattering (DIS) and Drell-Yan processes 
have been re-examined in the context of the gauge-invariant dehnitions 3, 3 3 • 

In semi-inclusive DIS and Drell-Yan processes, the quark TMD distributions produce 
the dominant contribution, while the contribution from the gluon distributions is power- 
suppressed Q. However, as a fundamental observable of the nucleon, the gluon TMD dis¬ 
tributions contain unique non-perturbative information, and contribute to a distinct class 
of semi-inclusive pro cesses [l^. Il^ liol 1^ . For example, a jet-jet correlation asymmetry in 
single-transversely-polarized-hadron-hadron scattering may reveal the role of the so-called 
gluon “Sivers function” Hence it is essential to have a thorough theoretical understand¬ 
ing of the gluon TMD distributions and related factorization theorems for gluon-related 
semi-inclusive processes. In general, however, QCD factorization for generic hard processes 
(e.g., heavy-flavor and jet productionj at hadron colliders is far more complicated than that 
for DIS and Drell-Yan processes [l3,|23, As a hrst step, we consider here scalar-particle 
production at hadron colliders through gluon fusion. We assume that the scalar particle 
(e.g. Higgs boson) couples to gluons at the leading order of an effective theory. All-order 
QCD corrections are taken into account in a factorization theorem. This process is similar to 
Drell-Yan production in that the observed particles have no hnal-state strong interactions, 
and thus QCD factorization can be studied in a similar manner. In this paper, we focus on 
the region of low transverse momentum where the TMD parton distributions are relevant. 

It is well-known that for a transverse momentum Aqqj-, Pf -C where being some 
hard scale in the problem (the scalar particle mass square in the present case), there exist 

large logarithms of type a™ln^”^~^ in high-order perturbative calculations. To have 

reliable predictions, we must re-sum these large logarithms 0, El S113. In the present 
framework, resummation can be achieved by solving the Collins-Soper evolution equation 
for the TMD parton distributions P, P . We will illustrate how to achieve the resummation 
for scalar-particle production at the relevant transverse momentum. 

One interesting example of the scalar particle is the standard model Higgs boson, whose 
discovery is among the most important endeavors for the future high-energy collider experi¬ 
ments. In the standard model, Higgs boson production at hadron colliders is dominated by 
the gluon fusion process: It couples to two gluons through quark loops, with the dominating 
contribution from the top quark loop. This coupling can be described by an effective theory 
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|28L I 29 L Isol. |3l| This effective theory is valid at the heavy quark limit, and in practice it 
works very well at the Higgs mass range of Mh < 2Mt jS^. In the past few years, there 
have been extensive discnssions abont Higgs prodnction in this effective theory approach, 
with the next-to-next-to-leading order corrections being calcnlated . Resnmmation of 

larg e log^ithm at low transverse momentnm was stndied by varions anthors in the literatnre 
3^ 37, sl, 3^, i3,|4l[|4^. In onr stndy, we nse the same effective theory for the conpling 


between the scalar and glnons. We start with the factorization analysis of scalar-particle 
prodnction at small transverse momentnm. For Higgs boson prodnction at large transverse 
momentnm, we conmare onr resnlt to the previons calcnlations with the resnmmation effects 


at one-loop order j35l . l3fiL 1371 l38l| . In onr analysis, large logarithms will appear in the TMD 
parton distribntions, and they are resnmmed by solving the corresponding Collins-Soper 
evolntion eqnation. The resnmmation coefficient fnnctions can be calcnlated directly from 
the factorization formnla. 

This paper is organized as follows. In Sec.H, we discuss TMD gluon distributions. First, 
we provide a gauge-invariant dehnition of the unpolarized distribution. We then calculate 
it to one-loop order in Feynman gauge. We show that the one-loop result is free of the soft 
divergence, and obeys the Collins-Soper evolution equation. In Sec. HI, we discuss in detail 
the connection between TMD and integrated parton (quark and gluon) distributions. In 
Sec. IV, we consider semi-inclusive scalar particle production at hadron colliders through 
the gluon fusion process, where an effective gluon-scalar vertex is introduced. We calculate 
the differential cross section and illustrate factorization at one-loop order. The relevant hard 
and soft factors are also calculated. In Sec. V, resnmmation of large logarithms is performed 
by solving the relevant Collins-Soper equations. We conclude and discuss future prospects 
in Sec. VI. 


II. TRANSVERSE-MOMENTUM-DEPENDENT GLUON DISTRIBUTIONS 


In this section, we study transverse-momentum-dependent gluon distribution, which is 
an important ingredient for scalar-particle production at hadron colliders. The discussion 
follows closely our previous work on the quark TMD distributions in semi-inclusive DIS 
[ 1 ^ . We hrst provide a gauge-invariant dehnition of the spin-independent TMD gluon 
distribution. We then calculate it in perturbation theory to one-loop order, and show that 
it obeys the Collins-Soper evolution equation. 

We focus the discussion on the unpolarized gluon TMD distributions. However, the 
results can easily be generalized to the polarized TMD gluon distributions 
with which one can study the various polarization asymmetries at hadron colliders. We 
leave those studies for future publications. 
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A. Gauge-Invariant Definition 


According to the previous studies QiQQlIi, a gauge invariant, spin-independent 
TMD gluon distributions can be defined through the following matrix element 


xg{x,k±,n,xC,p) = 


P+(27r)3 

'^i^;+„(r.n)£L(r.a)As=(o,ox)P'‘+('))i^’ 


X ^-^77-^> (1) 

where we have included a soft factor following the TMD quark distributions [l^ . In Sec. 
IIB, we will discuss the definition of the soft factor and its contribution. In the above 
equation is the gluon field strength tensor, — d’^A^ — gfabcA^A’^. We use 

the convention for the gauge coupling + igA^. Light-cone components are defined 

as = {k^ ± k^) /a/ 2. Therefore F"*" is the light-cone momentum of the hadron, and x is 
the momentum fraction carried by the gluon, while k± is the transverse momentum. £„ is 
the gauge link from the past. 


( 2 ) 


C)o) = Fexp J dXv ■ A{Xv + ^) 

xPexp^-igj dr]± ■ A_i{-oo,g_i_)^ , 


where A^ = A^P is the gluon potential in the adjoint representation, with = —ifabc- 
Four-vector v is an off-light-cone vector v = {v~,v^,v± = 0) where v~ S> n’*'. In principle, 
the gauge link is best taken along the light-cone direction. Such a light-cone gauge link, 
however, introduces the light-cone singularities The g aug e link is therefore chosen 

to be slightly off the light cone to regulate these singularities j^. With a non-light-like v, 
the TMD distribution depends on a new scalar = (2v ■ F)^/n^. In light-cone coordinates, 

= 4:{P~^)‘^v~/v~^ and is proportional to the hadron energy. Energy evolution for the TMD 
distributions is governed by the so-called Collins-Soper equation [l|. /r is the standard renor¬ 
malization scale. We use the dimensional regularization and modified minimum subtraction 
(MS) to take care of ultra-violate divergences. 

The gauge link above can be derived in the same way as that in the quark distributions |l^ . 
Its direction is chosen to start from —oo, indicating its origin from initial-state interactions 
present in hadron-hadron scattering. This is similar to the gauge link in the TMD quark 
distributions appropriate for Drell-Yan production l3|. [One can, of course, choose a gauge 
link pointing to the future to define a similar distribution relevant, for example, for DIS.] 
In nonsingular gauges, the transverse gauge link vanishes because gauge potentials fall off 
rapidly at space-time infinity [l^ . In a singular gauge, e.g. the light-cone gauge, gauge 
potentials do not vanish at space-time infinity, so we have to include the transverse gauge 
link to guarantee gauge invariance Q. In the following calculations, for simplicity, we will 
work in Feynman gauge in which the transverse gauge link does not contribute. 

We are interested in only the leading-twist contribution. Therefore, we keep only the 
leading power in 1/C^, and take the ^ oo limit whenever it is convenient and there are 
no singularities. 

Feynman rule for the vertex connecting the probing gluon and the gauge link is 
shown in Fig. 1. The first term comes from the gauge link, while the second one from the 
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FIG. 1: The Feynman rule for the TMD vertices. The first comes from the eikonal line and the 
second comes from the non-ahelian term in the field strength tensor. In this calculation, the sum 
of the two will be lumped together and labelled by the gauge link. 





FIG. 2: Virtual corrections to the TMD gluon distribution at one-loop order. The mirror diagrams 
are not shown. 


nonlinear part of the tensor. In the following calculations, we present the results for the 
contributions from these two vertices as one term and label it as the gauge link contribution. 


B. One-loop Calculations Without Soft Factor 

In this subsection, we calculate the unpolarized TMD gluon distribution in a gluon target 
in perturbative QCD to one-loop order, ignoring the soft contribution. According to the 
definition, the distribution is normalized to 

g{x,k±) = 6{x - l)5‘^{k±) , (3) 

at leading order. At one-loop order, we have both virtual and real corrections. We regulate 
the infrared (soft) divergences using dimensional regularization. Moreover, we introduce the 
off-shellness for the initial gluon, < 0, to regulate possible collinear divergences. 

The one-loop virtual diagrams are shown in Fig. 2, which are the self-energy diagrams 
for the gluon wave function and from the gauge link, and the vertex corrections. The 
contribution from Figs. 2(a), (b) and (c) has the form 

g{x, /c±)|fig.2(a-c) = S{x - l)5‘^{k±){ZG - 1) , (4) 
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FIG. 3; Same as Fig. 2 but for the real contributions. 


where 


Zq — I + 


47r 



31 


10 


+ - -N., 


(5) 


A factor of 2 has been included to account for the contributions from the conjugate diagrams. 
Here Ca = and Nf is the number of the active quark flavors. Minimal subtraction has 
been made to remove the UV divergence. Collinear divergence is indicated by the dependence 
on the off-shellness of the gluon. 

The contribution from Fig. 2(d) has the similar form, 6{x — l)5‘^{k±){Zw — 1), with 


Zw = 1 + 


ttgC/ 

2ti 


'2 1 

-h 7E + In — 

Gr 


( 6 ) 


where 1 /cir (eiR = 4 — d) pole comes from the soft divergence. 

The contribution from Fig. 2(e) is 5{x — l)S^{k±){Zv — 1) with 


Zy — 1 — 


C^sC/ 

271 
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dvr/i^ 


2hF 
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dvr/i^ 


12 


(7) 


Here we have taken the limit = {2v ■ P)^/n^ is much larger than any soft scale (e.g. —p^). 
The result shows explicitly the double logarithmic dependence of the TMD distribution on 
These double logarithms can be resummed using the Collins-Soper evolution equation, 
which will be discussed in Sec. HE. The total contribution from the virtual diagrams is 


5((a;, fc_L,h,0 [virtual = ^(a;-l)d^(fc_L)(ZG + Zvy + Zy -3) . (8) 


Again the virtual correction contains soft divergences. 

Now we turn to the real corrections shown in Fig. 3. The contribution from Fig. 3(a) is 

/ 7 \ I _ A 1 

fc±)|fig.3(a) - ^2 kl-x{l-x)p'^ 

Fig. 3(b) contributes 

9(r. fcr. xC) k,.3(., = ^ 

( 10 ) 


X 


X 


(1 -x). 


, 1 h' , C r/ .X 

Zd n, \ rv \ 


1 — X . X 

-^ + x{l-x) + - 


(9) 
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FIG. 4: Soft contributions to the TMD gluon distribution at one-loop order. 


where a factor of 2 has been included to account for the mirror diagrams. The plus function 
follows the dehnition in 3|- The limit of > cxo has been taken whenever there is no 
divergence. The contribution from Fig. 3(c) can be calculated similarly, 


k±) |fig.3(c) 


S' 

27r2 k‘j_ 




1 ) 


( 11 ) 


The above contributions contain soft divergences as indicated by the l/k^ term as /c^ ^ 0. 
We have to keep the pre-factors of these terms up to 0{e), which can lead to hnite contri¬ 
butions after Fourier transforming into the impact parameter 6-space. The soft divergences 
from the real and virtual diagrams will eventually cancel, as we shall see. 


C. Soft Factor 


As in the case of TMD quark distributions, there are soft contributions from Figs. 2 and 
3, signaled by soft divergences. These contributions can be extracted using the well-known 
Grammer-Yennie approximation jd^. Here we follow the same procedure as for the quark 
distributions 113 , and dehne the soft factor 

S{b,s,p) = X) -oo)Cl^,,^{-oo] 6±)£^,ab(0; -oo)Cybc{-oo] 0)|) , (12) 

where v is another off-light-cone vector, v = (h“,h’'',0) with ^ and p is dehned as 
= v~v^/v^v~ . The same soft factor will appear in the factorization of semi-inclusive 
scalar particle production we shall see in Sec. IV. 

At one-loop order, we have diagrams shown in Fig. 4 which are the same as those 
considered in except for the color factor. Here we use dimensional regularization for 
the infrared divergence instead of a gluon mass. The hnal result is the same because the 
soft factor is free of infrared divergences. 

The contribution from Fig. 4(a) is 

Asoft5'(^, /c_L)|fig.4(a) = -25{x - l)5‘^{ki_){Zw - 1) , (13) 

where Zy/ is in Eq. ©. The contribution from Fig. 4(b) is 

Asoft5'(a;, k±) |fig.4(b) = 6{x- 1)(5^(A;_l)I n f^ -h 7 e -h In ■ ^ 


27r 


Hr 


dvr 


(14) 


where the UV divergence has been subtracted using MS scheme. The contribution from 
Fig. 4(c) is 

Asoftfi'(a:, A;_L)|fig.4(c) = b{x - 1)^^^-^^ ■ (15) 


TT^ 
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Fig. 4(d) gives, 


(16) 


As„ri9(i, *:i)|iig,4(d) = -i(x - ■ 

The individual diagrams all have infrared divergences. However, they all cancel eventually. 
For example, the IR divergence from Fig. 4(a) cancels that from (c), while that in (b) cancels 
that in (d). The cancellation can best be seen from the impact parameter 6-space expression. 
After Fourier transformation, we get 

Asoftg(x, b, fi, x(, p) = S(x - In [lnp^-2] (17) 

We thus conhrm that it is the same soft factor as in the quark distributions Q , except for 
the color factor. 

We note that there are different dehnitions of the soft factor in the literature [ilElil. 
The difference can be viewed as the scheme dependence. Final physical results do not depend 
on a particular dehnition. 


D. TMD Gluon Distribution at One-Loop 


Summarizing the results from the last two subsections, we get the gluon TMD distribution 
at one-loop order, 

g(x,k±,p,x(,p) = S(x - l)S^(k±) [Zq + Zy - 2 + {Zw - l)iMp^ - 1)] 


+6{x - 1) 


p 


27r2 kl 
O^sCa 1 


In - In -h 1 

k. 


X 


H-h x(l — X) 

X 


(18) 


7 r2 k^ — x{l — x)p‘^ [(1 —a;)+ 

The individual terms have soft divergences which cancel out in the sum. In the impact 
parameter 6-space, we dehne 


g(x,b,ii,x() = I d‘kxe'^^*^g{x,kug,xC) 


(19) 


After Fourier transformation. 


g{x,h,p,xC) = 


a.gC/ 


vr 


In ■ 


-p'^b" 


- 27b 


X 


+ 


Oi.qC/ 


TT 


X 


In 


1 — X x{l — x) 


(1 — x)+ X 

1 \ , 1 

-|- In 


^ - X , 

H-h x(l — x) 


x{l — x) 


X 


X 


+ x(l — x) 


+ 5{x - 1) {Zg - 1) + 


UsCa 

2tt 


--In^ 

2 


(^62 


,27b 


In 


—p%‘^ 


plE 


In 


c,v 


+ In (in/ - 1) + 1 hX - f! - 

4 ^ ^ 2 p2 2 


4 

^2 n' 
2 2 


( 20 ) 


The 1/eiR, poles have now disappeared. However, there are still collinear divergences as 
indicated by the dependence on ln(—p^). More interestingly, we have double logarithmic 
dependence on <^2, the energy of the parenting gluon. We study this dependence in the next 
subsection. 































E. Collins-Soper Evolution in Hadron Energy 


From the above result, we see that the TMD gluon distribution depends on the energy 
of the parenting hadrons, through the variable = (2nP)^/n^ ~ 2(P+)^n“/n+. The en¬ 
ergy evolution of the TMD parton distribution is controlled by the Collins-Soper evolution 
equation [ij. In impact-parameter space, the Collins-Soper equation reads as 


d 

b, xC, ft, p) = (K, + Gg)g(x, b, xC, p., p) , 


( 21 ) 


where K and G are soft and hard evolution kernels, respectively. 

It is easy to check that the one-loop TMD gluon distribution indeed satisfies the Collins- 
Soper equation. The sum oi K + G can be extracted. 


Kg + Gg 


asGA , 27B-3 

-In-:—e 2 

TT 4 


( 22 ) 


The soft part K can be calculated from a definition similar to the quark case 


Hii, 


= 




In ^:—e" 
TT 4 


(23) 


and thus the hard part Gg can be solved from the sum. The K and G obey the renormal¬ 
ization group equation. 


dKg dGg a,GA 

= -iKg = -2 - 

on on TT 


(24) 


These evolution equations will be used to resum the large logarithms in the cross section in 
Sec. V. 


III. TMD PARTON DISTRIBUTIONS AT LARGE TRANSVERSE MOMENTUM 


As emphasized in HI. it is nontrivial to generate the integrated (Feynman) parton dis¬ 
tributions by integrating out the transverse momentum in the TMD distributions. The 
main difficulty is that new ultraviolet divergences emerge from the transverse momentum 
integration. However, there is another way to connect the two distributions j^: When the 
transverse momentum becomes large or the impact parameter becomes small, part of the 
TMD distributions are calculable from perturbative QCD, because at least one hard gluon 
is needed to generate the large momentum. The hard-gluon radiation leads to power-like 
behavior, e.g. for the unpolarized TMD distributions. Here Feynman parton distri¬ 

butions enter as the non-perturbative input in the QCD factorization analysis. 

We consider the Fourier-transformed version of TMD distributions at small b. The fac¬ 
torization theorem reads Hi, 


fi{x,b,n,xC,p) = Y^ !—Gi 
j ^ 


i/j 






(25) 


where i, j denote the flavor of the partons—quarks and gluons. The integrated distributions 
fj on the right-hand side depend on the longitudinal momentum fraction x and the scale Jr. 
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Cijj are the coefficient functions and are calculable from perturbative QCD. In the following, 
we will verify the above formula and extract the coefficient functions at one-loop order. 

We choose an on-shell quark target (with mass m) and an off-shell gluon target (with 
off-shellnes p^) to regulate the collinear divergences. We need to calculate both the TMD 
and integrated parton distributions up to one-loop order. The collinear divergences in both 
distributions must be matched, and the subtracted coefficient functions are free of infrared 
divergences. 

The coefficient functions Ci/j have perturbation expansions in terms of the strong coupling 
constant At leading order ( 0 ;°), it is easy to see. 


C 


( 0 ) 

9/9 

( 0 ) 

g/q 


C 

C 


( 0 ) 

9/9 

( 0 ) 

9/9 


6{x - 1) , 

0 . 


(26) 


At one-loop, all receive non-trivial contributions. Cg/q has been calculated before in Q, 
where a gluon mass was used to regulate the soft divergences. We have checked this calcu¬ 
lation by using dimensional regularization and obtained the same result. For completeness, 
we quote the result here. 


Cq/q (a;, p) 

1 + 


Qfc 




1 — X 






+6{x - 1) 


T-lnp')ln^e - In 




,27b-1 


2) + Tl^ 


(27) 


To calculate we need the integrated gluon distribution at one loop. 




'gig 

A 


IT 




1 — a; 


X 


-f a;(l — x) ) In 


— In —-- 

X x(l — X) 
1 


. 31 5Nf 


x(l — x) j ’ 


(28) 


where Vgg{x) is the gluon splitting function. 


'Pggi^) 


X 


(1 -x). 


1 — X 


X 


-F x(l - x) + S(x - l)/3o , 


(29) 


with /3o = (11 — 2/3Nf)/12. From the above and one-loop TMD gluon distribution, the 
coefficient function for the gluon-gluon term is 


*^9/9 p) 


a.,C' 


A 

7T 


Vgg{x) In -r^e -F 5(x - 1) 


4 4 


h^jjL 


^0 + lnp - ^ 1 


3 27b 


.27b 


71 ^ + 7 


( 30 ) 


where the collinear divergences of the form ln(—p^) have been cancelled. 
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FIG. 5: Quark splitting to gluon contribution to Cgjq (a) and gluon splitting to quark Cq/g (b). 


The Feynman diagrams for the qnark-glnon and gluon-quark splitting coefficient functions 
Cg/q and Cq/g are shown in Fig. 5(a) and (b). The TMD gluon distribution in a quark target 

is 


k±) |fig. 5 (a) 


ttsCp 

27r^ 


1 -I- (1 — xY 1 


2 a;(l — x)rri^ 


X k\ + m? {k\ + x‘^m?y 


After Fourier-transforming into the impact parameter 6-space, we have, 

OisCp 


(a;, 6) I fig,5(a) 


27r 


l + i„ ( 4 _ 2(1-:.) 

\mVx^ 


X 


X 


The corresponding integrated gluon distribution is, 

UgCp 


g{Xi Y) |fig. 5 (a) 


27r 


1 -1-(1 — x)^, 2(1 — x) 

-In —- X 


X 


m?x‘^ 


X 


Combining the above, we find the coefficient function. 




2i 


1 -I- (1 — xY , / 4 

In I \ + X 


X 




Similarly, we hnd the coefficient function for hnding a quark in a gluon target. 




(y.g 

47r 


+ (1 - xf ) 111 ( ) + 2:.(1 - :.) 


(31) 


(32) 


(33) 


(34) 


(35) 


These calculations can certainly be carried out to higher orders in ag 


IV. SCALAR-PARTICLE PRODUCTION THROUGH GLUON FUSION 

It is difficult to study the gluon TMD distributions in semi-inclusive DIS because gluons 
do not have direct couplings with photons, and their contributions are power suppressed 
[T^ . However, we can access these distributions at hadron colliders through, for example, 
the heavy-quark pair or di-jet production. In these processes, the total transverse momentum 
of the pair or the di-jet provides information on the TMD parton distributions in the initial 
hadrons, and the gluon distributions could dominate at some kinematics. 

An important question we have to address is the factorization of the QCD radiative 
corrections. Only when the factorization is established can we safely extract the TMD parton 
distributions from data. In this section, we will consider QCD factorization for the gluon 
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FIG. 6: Vertices for the scalar particle coupling to the gluons. All labelled momenta are in. 


initiated processes, with scalar-particle production as an example. The scalar particle here 
has a direct coupling to the gluons through an effective vertex. The transverse-momentum 
distribution of the produced scalars reflect partly the transverse-momentum dependence of 
the gluon densities in the nucleon. Our study can be extended to other gluon initiated 
semi-inclusive processes at hadron colliders. 

The effective lagrangian for the gluon-scalar coupling is, 

, ( 36 ) 

where <h is the scalar field and gcj, is the effective coupling. For the standard model Higgs 
boson, the effective coupling can be derived from the full theory [ 2 ^ |2^ . The above effective 
lagrangian has also been used to study the gluon saturation in nuclei. 

From the above lagrangian, we can write down the basic vertices for the scalar particle 
coupling to the gluon potentials, shown in Fig. 6. The Feynman rule for the coupling with 
two gluons in Fig. 6(a) is 

'^9(j)^ab (^1 ■ . ( 37 ) 

And the coupling with three gluons in Fig. 6(b) reads 


9<j>9sfabc [(^1 ^2}p9uF ff (^2 k'i)p.9pv T (^3 ^l)u9pp\ ■ 


( 38 ) 


Finally, the coupling with four gluons shown in Fig. 6(c) is 


^9(l>9s { fabefcde (.9p,p9 u'y 9p,'y9yp) 

T facefbde {,9^139p-f 9p,'y9i'p') 

T fadefbee (.9P.1'9pr/ 9pp9i''y)} 


( 39 ) 


Production of a single scalar particle has a cross section linear in g"^. We are interested in 
higher-order QCD corrections in as- 

Scalar-particle production at low transverse momentum depends on the TMD gluon dis¬ 
tributions of the incident hadrons. A new feature for semi-inclusive processes is the depen¬ 
dence on the soft factor 0 , 0 . In a factorized form, the cross section can be written as 
a product of the TMD parton distributions (or fragmentation functions), the soft factor S, 
and the hard factor H [l|, [ 13 , 


dMM^P±,y) 

(PP±dy 


(To J d‘^ki_Ld‘^k22_d‘^i±xig{xi,ki±,xiCi, Id, p)x2g{x2,k2±,X2C2, y, p) 

xS{i^,p,p)H{M\p,p)6\h^ + k2± + e^-P^) , ( 40 ) 
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( 41 ) 


where cxo is the leading-order scalar-particle production from two gluons, 

_1 

“ 64 1 - e/2 ■ 

Here the factor (1 — e/2) in the denominator comes from the polarization average of the 
initial gluons. The hard scale M here denotes the mass of the scalar particle, and y and P± 
are its rapidity and transverse momentum, respectively. At low-transverse momentum, the 
longitudinal-momentum fractions Xi and X 2 for the two incident gluons can be approximately 
related to the scalar particle’s rapidity y through xi = \/M‘^ jSe^ and X 2 = \/ M^/ 
where S is the total center-of-mass energy squared S = (Pi -|- P 2 )^. Ci ^md (2 are dehned 
as Cl = 4(n ■ Pi)/v‘^ and Cl = 4(h • P 2 )/h^, and v and v are dehned in the last section, p as 
dehned before, is a parameter to separate gluon contributions to the soft and hard factors. 

The above factorization result is accurate at the leading power in P^/M^ at low transverse 
momentum. In the following, we will examine this factorization formula up to one-loop order 
using the perturbative TMD gluon distributions from the previous section. Through explicit 
computation, we obtain the one-loop result for the hard factor H. For arguments toward a 
factorization to all orders, we follow the discussion in [iL ll^. 

In principle, the above factorization formula works only at the low transverse momentum 
region, and it breaks down at P± ~ M, where one can no longer neglect the power corrections 
of P]_/M‘^. However, it is useful to extrapolate the above factorization to all Pi, and convert 
it to the impact-parameter space. After Fourier transformation, the diherential cross section 
can be written as. 


dMM^Px,y) 

d?Pi_dy 

where we dehne 


o-Q 


W? 




W{x,,X2,b,M^) , 


(42) 


W{xi,X2,b,Q^) = xig{xi, b, xiCi, p, p)x2g{x2,b, X 2 C 2 , p)S{b, p, p)H{M^, p, p) . (43) 

The convolutions in the transverse-momentum space now reduce to products in the impact 
parameter 6-space. 


A. Factorization At One Loop 


In this subsection, we calculate the scalar-particle production to one-loop order to verify 
the factorization, and to extract the hard factor H. The gluon TMD distributions and 
soft factor have been calculated in Sec.H, where the gluon target has been put off-shell (by 
nonzero —p^). In this section, it is more convenient to demonstrate the factorization by 
using the on-shell gluon target (p^ = 0). It is straightforward to extend the results in Sec.H 
to the on-shell case. In the new scheme, the perturbative TMD gluon distribution is. 


g{x,b, xC,p,p) 


5{x — 1) 

IT 


'Pggi^) 


+6{x - 1) 


3 , C' 7r2 + 7' 
■4 'V —— 


\np + po - 


In 


Anp^h^ 

4 4 



( 44 ) 
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FIG. 8: Real diagrams contribution to the scalar particle production through gluon fusion: part I. 


Here, l/ejR indicates collinear divergence. The soft factor S{h,pi,p) was defined in Eq. (fT^ . 
At one loop, the result in impact parameter space is 

(2-lnp2) . (45) 

These results will be used to verify the above factorization formula. 

At leading order, the Born diagram for the production cross section is shown in Fig. 7. 
The calculation is straightforward, and we get 


This leads to 

h) = 5{xi - 1)5 {x2 - 1) , (47) 

which can be clearly reproduced by the leading-order gluon distributions in gluon targets 
At next-to-leading order, we have both real and virtual contributions. The real diagrams 
are shown in Figs. (8-11). Figure 8 represents the contributions from the gluon radiation from 
either gluons, while Fig. 9 for the interference between these radiations. Figure 10 includes 
the diagrams involving the three-gluon coupling with the scalar particle. The contribution 
from Fig. 8 is 

5(a;2 - 1) + (xi ^ 0 : 2 ) 1 (48) 


da 


d?P^dy 


= (Tq- 


fig.8 


Cac^s h " 


TT^ 


Xi 


1 -a;i xi 

-h xi(l -xi) + — 

Xi 2 


In the above result, we have kept the pre-factor up to 0{e), because 1/Hf will lead to infrared 
divergences when Fourier-transforming to the impact parameter space. The contribution 
from Fig. 9 is. 


da 


d'^P^dy 


= CTq- 


^A^s R 


fig.9 


vr^ 


Pi 


Xi 


(1 -Xi). 


h 

2 


1 

- 1) 


5{x2 - 1) 


+ (Xi ^ X 2 )} 


(49) 
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(a) {b) 

FIG. 9: Real diagrams: part II. 






FIG. 10: Real diagrams: part III. 


and the contribution from all diagrams in Fig. 10 is, 


da 


d^P^dy 


= 0 ' 0 - 


fig. 10 


CAOis 

2 p2 


71 


xf + Xi 


5{x2 - 1) + {xi ^ X 2 ) 


(50) 


The diagrams in Fig. 11 are power-suppressed at low-transverse momentum, although their 
contributions are leading when the transverse momentum is on the order of M. The sum of 
the contributions from all the real diagrams is 


da 


d^P^dy 


= (To 


real 


7r2 


Xi 


Xx 


H-h Xi{l - Xi) 


1 

._ln_5(a;,-l) 


(l-a;i)+ xi 
5{x2 - 1) + ^ X 2 ) 


(51) 


The above agrees with a corresponding result in Ref. j37[ . After Fourier-transforming to the 
impact parameter space, one hnds. 




(XsCa 


TT 


-7 e + In 


<;IR 




5 { x 2 — l)xi 


Xl 


(1 -a^i) + 


-Fxi(l - Xi + ^——) ) -h (xi ^ X 2 ) 


Xl 


-h(5(xi - 1 ) 6 { x 2 - 1) 


2 m2 \ 1, 2 / m2 

In 1 ' 1 


. /I 2 ^ 

LGr Gr V 47r/i"= / 2 


47rp2 


— In2 




M2 


Te 


+7pln^^ + ^- 
' 47r/T2 ^2 12 


( 52 ) 


The double and single poles above correspond to infrared divergences only. 
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FIG. 12: Virtual diagrams contribution to the scalar particle produetion through gluon fusion. 


Now we turn to the virtual diagrams shown in Fig. 12. These diagrams contribute to the 


cross section as 


33,0 


asCA 


W{xi,b,M^)\Artne.\ = 5(a;i - l)(5(a;2 - 1) 

TT 

„2 ^^2 


4 2 /, m2 \ 1 , 2 / m2 

In 1 ' m2 


-2 

e 




, M2 Ye 77r^ 


47ryu2^ 

(63) 


The l/e2 pole above cancels that from the real diagrams. The above result has ultra-violet 
divergences which comes from the composite gluon operator in the effective coupling. This 
divergence can be cancelled by a charge renormalization contribution, 


IT(a:„6,M2)|,harge = ^2/3o ^ " 1) • (54) 

where the extra factor is introduced so that renormalization scale is set at j^jy = M2. The 
gluon wave function renormalization constant Za = 1 to all orders in perturbation theory in 
the present scheme. By summing up all the contributions, we get W{xi,b, at one-loop 
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order, 


O^sCa 

71 


XlVgg(Xl)6(X2 - 1 ) 


+ 6{xi - 1)6 {x2 - 1 ) 


2 /?oln 



- 7i? + In 


ATTjJ?})^ ) 


+ (Xl 




( 66 ) 


The soft divergences have been cancelled between real and virtual diagrams, and there are 
remaining collinear divergences. 

Now we can verify the factorization formula in Eqs. dm), dini), using the one-loop results 
in Eqs. dH, dlSI)- After the subtractions, we obtain the hard factor at one loop. 




Q^sC'a 


TT 




3 2 1 1 2 2 2 

- In p + - In^ p-^ + 7r^ + - 
4^ 8 ^ 2 


( 66 ) 


where a special coordinate system has been chosen: xf(i = x^C^ = pM"^. Not only the 
collinear divergences in the cross section have been cancelled by the gluon distributions, 
the 6-dependence has also been entirely isolated in the parton distributions and soft fac¬ 
tor. The hard factor H depends only on the hard scale M^, the factorization scale p, and 
the parameter p. Thus, the factorization conjectured in Eq. (HOD is verihed to one-loop 
order. Furthermore, according to the general arguments provided in m, all-order QCD 
corrections to scalar-particle production can be factorized in a similar way. 


B. Factorization at Small b <C 1/Aqcd 


When b is small, there is the standard QCD factorization in terms of Feynman parton 
distributions for W (xj, b, M^) [^. In this subsection, we will try to recover this factorization, 
using the results in the last subsection. 

According to [^, the cross section IE(xj, 6, M^) at small b can be written as. 


d^i d^2 


Xi 


,X2 


W{xi,b,M ) = / -— p)^ 29{^2, p)Cgg{ — ,b, p, M )Cgg{ — ,b,p,M ) , (57) 


6 6 




6 


where g{x, p) is the usual gluon distribution, and the coefficient functions Cgg is perturbative. 
They both depend on the factorization scale p, but the dependence cancels out in the 
product. 

Using W{xi,b, M"^) from the last subsection, we can calculate the Cgg to one-loop. The 
perturbative gluon distribution up to one-loop is 

g{x) = 6{x-l) + '^^Vgg{x)(^-‘^ + -fE + \n^^ . (58) 


Subtracting the above from the one-loop W{xi, b, M^), we have 


Cgg{x,b,p,M^) = 5(a;-l)-F 


+6{x - 1) 


a.^C/ 


TT 


XVgg{x) In 


Po\n 


b‘^p‘^ 
1 


3 -27b 


fb^M^ 2-v \ 2 




,27b 


+ T 
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The /i-dependence of the Cgg function is controlled by the same evolution equation as for the 
gluon distribution, with the opposite sign. When is large, one needs to re-sum the large 
double logarithms to all orders in perturbation theory. This is difficult to do systematically 
in the present formalism. However, this can be done straightforwardly in terms of the TMD 
factorization in the previous subsection. 


V. RESUMMATION OF LARGE DOUBLE LOGARITHMS 

The one-loop result for scalar-particle production shows that there exists large double 
logarithms in the hard scale (e.g., see Eq. (j^). These large logarithms appear at every 
order of the perturbative expansion, a^ln^^M^/P^ where m <n and P± is the transverse 
momentum. In order to make reliable predictions, one has to re-sum these large logarithms. 
In this section, we will follow the Collins-Soper-Sterman (CSS) framework jsf to perform 
the resummation. The results can be used as a double check of the resummation studies in 
the literature for the standard model Higgs production at hadron colliders 


A. Collins-Soper-Sterman Resummation at Small b 


According to the QCD factorization, the Fourier transformation of the cross section for 
scalar-particle production is 

W{xi, b, M^) = xig{xi, b, p, pM^, p)x2g{x2, b, p, pM'^, p)S{b, p, p)H{M^, p, p) (60) 

In the above equation we have set x^^f = X 2 C 2 = in a special frame of coordinates. 
The dependence of W can be studied through the following differential equation 


dW{xi,b, 

9 In 


= {K + G')W{xi,b,M^) , 


(61) 


where K and G' are soft and hard evolution kernels. The soft part K depends on the scale 
1/6^ and the renormalization scale p, while G' depends on the hard scale and p. K can 
be calculated from the Collins-Soper evolution equation for the TMD gluon distributions 
discussed in Sec.HD. G' contains the contributions from the gluon distribution as well as 
from the hard factor H. From the result of the previous sections, the sum of K and G' at 
one-loop order is. 


K + G' = 


UsGa, 


In 


TT 


pi'yE—'iPo 


(62) 


where the p dependence between various terms cancels out. Our previous result for K at 
one loop is 


K = Kg = 


In , (63) 

71 4 

which is valid when b is small enough. At large b, the perturbative expansion breaks down. 
The hard part G' is. 


G' = 


OLsGa 


TT 


1 ^ 

In — - 2/3o 

/i2 


(64) 
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Both the soft and hard parts K and G' also obey the renormalization group equation j^, 

dK dC 


9 In /i d In fi 


= -iKg • 


( 65 ) 


Resummation is done here by solving the above equations. 
Integrating over InM^ and In/i^, one hnds 


where the Sudakov form factor is 

c _ 

<^Sud — —5 

Jcf/b^ h 


inl^ 

/i 


fi) + B{Ci, C21 /i) 


( 66 ) 


(67) 


Here C\ and are two parameters of order one. The functions A and B can be expanded 

00 . 00 . 

perturbatively ag, A = Hh) (^)* and B = Y (^)*- They are related to the Collins- 

i=l i=l 

Soper evolution parameters by the following equation, 


MC,,^i) = \iKM + \liAK{Ci,g(tA) , 
B{CuC 2 ,^) = -K{Cug(p))-G'(l/C,,g(p)) . 
From Eqs. ()6dl64j] . we get the hrst expansion of A and B functions, 

= Ca, = Ca 


In / 


4C| 


-2/?o 


( 68 ) 


(69) 


To complete the resummation, we use the factorization of W{xi,h,Cl/Cl/h^) in terms of 
the integrated parton distributions discussed in Sec. IVB, 


W{xi,b, 


Cl 

Clb^ 


d^i 




]C, 


iC'civ 




i=0 


,6’ 'Clb^' 

x6^(6,h)6^(6,h) , (70) 

where C functions also have perturbation expansion in terms of C = Y (^) • From 
the results in Sec. IVB, we get the hrst two terms of Cgg, 

d? = Hx - 1) 

Gs = lm(x) In f 


-F5(x - 1 ) 




4C|' 


3275 


TT 


+ T 


(71) 


It is easy to check that the above C function for the gluon-gluon term can also be calculated 
from the following equation 0 






^*39(^5 ^2^2 ’ ^ ^ xCggix^b ,/ig, p 1 P Q2jg2'’P j ^ \ I Pgi P)dd { q2Jj2’’P9 '’’ ( 72 ) 
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with Cgg from Eq. dsn}, and S and H from Eqs. (03), dsni) respectively. Similarly, the contri- 
bntion from the qnark distribntion in the factorization form of Eq. dzni) can also be calcnlated 
nsing the above eqnation with Cg/g in Eq. Since = 0, it is straightforward to get 

the expansion of Cg/g np to one-loop order, 




0 

Cf 

—X 


ri + (i 



X 



+ X 


(73) 


Onr hnal resnmmation resnlt is the combination of Eqs. (66), (67), (70), (71) and (73). 


B. Comparison with Previous Calculations 

In the literatnre, the above coefficient fnnctions have also been calcnlated by varions 


anthors j35l. I36L 1371 l38| . They are normally extracted from the comparison between the 


hxed-order results and the expansion of the resnmmation formula, and expressed with the 
so-called canonical parameters, with Ci = 2e“'^® and C 2 = 1. Choosing these specihc values, 
our results for A, B, C at one-loop order are. 


AW = Ca, BW = -2/3„C^, = Ca^S(x - 1) . 


(74) 


AC) and i?C) agree with the previous calculations jSS], |3a |37| , while lacks an additional 
term 11/4. The difference comes from the effective coupling between the Higgs boson and 
gluons 113, im 


9h = - 


a.. 


3V TT 


1 + 


11 a. 




TT 


(75) 


where V is the vacuum expectation value of the Higgs held. The second term will contribute 
at one-loop order in ag and modify the result of ()53|1 . Taking this into account, our 


results coincide with the previous calculations |3a |37 . 

In , the full of dependence on Ci and C 2 of the coefficient functions were also calculated 
for the standard model Higgs boson production. Our result on i?d) agrees with theirs, while 
does not 


VI. CONCLUSION 

In this paper, we have studied scalar-particle production at one-loop order to examine 
the factorization of the gluon initiated semi-inclusive processes at hadron colliders. We in¬ 
troduced and calculated the gauge-invariant TMD gluon distribution at one-loop order, and 
established their connections with the integrated parton distributions at small h. We verihed 
the factorization theorem at one-loop order, where the cross section can be decomposed into 
products of TMD parton distributions, soft and hard factors. The large logarithms in the 
cross section were resummed following the CSS formalism, and the coefficient functions were 


^ After correcting some errors in |33, these two calculations agree on 
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obtained. The factorization to all orders in QCD conpling is assnmed following the general 
argnments in [l[ |l^ . 

The present stndy can be easily extended to the polarized glnon TMD distribntions and 
their contribntions to the spin asymmetries for the semi-inclusive processes. Scalar-particle 
production can also be generalized to other production processes, e.g., heavy-quark pairs 
and heavy quarkonium production, di-photon and di-hadron and photon-hadron production, 
and di-jet correlations. We will carry out further studies in future publications. 

The present study of factorization for semi-inclusive processes can be extended to two- 
loop order, which can be used to compare with the NNLO calculation of the standard model 
Higgs production. It will be interesting to calculate the resummation coefficient functions 
from the present approach, and compare with the previous results obtained by the expansion 
method 3^. 

As a hnal remark, we consider the small-x property of the gluon TMD distributions. 
From the results in Sec.II, we found that the TMD gluon distribution at small x has the 
similar behavior as the integrated gluon distribution, where the small x resummation is 
usually necessary. The so-called BFKL evolution equation is usually used to make the 
resummation. It will be interesting to examine whether the TMD gluon distributions obey 
such evolution equation at small x and to explore the connection between the present 
approach and the so-called fct-factorization approach [HJ. Phenomenologically, the small-x 
effects are important for the Higgs boson and vector-boson production at LHC IHT 
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